A minimal model for the inelastic mechanics of biopolymer networks and cells 
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Live cells have ambiguous mechanical properties. They were often described as either elastic 
solids or viscoelastic fluids and have recently been classified as soft glassy materials characterized by 
weak power-law rheology. Nonlinear rheological measurements have moreover revealed a pronounced 
inelastic response indicative of a competition between stiffening and softening. It is an intriguing 
question whether these observations can be explained from the material properties of much sim- 
pler in-vitro reconstituted networks of biopolymers that serve as reduced model systems for the 
cytoskeleton. Here, we explore the mechanism behind the inelastic response of cells and biopoly- 
mer networks, theoretically. Our analysis is based on the model of the inelastic glassy wormlike 
chain that accounts for the nonlinear polymer dynamics and transient crosslinking in biopolymer 
networks. It explains how inelastic and kinematic-hardening type behavior naturally emerges from 
the antagonistic mechanisms of viscoelastic stress-stiffening due to the polymers and inelastic flu- 
idization due to bond breaking. It also suggests a simple set of schematic constitutive equations 
which faithfully reproduce the rich inelastic phenomenology of biopolymer networks and cells. 

PACS numbers: 87.16.ad, 83.60.-a, 83.10.-y 



I. INTRODUCTION 

The mechanical properties of cells strongly influence 
processes in living organisms on all length scales, from 
the locomotion of single cells to the expansion of whole 
tissues during breathing. Therefore, studying the ma- 
terial properties of cells has a long tradition [1 . One 
of the major lessons learned so far is that cells are nei- 
ther solid nor fluid, but can tune their mechanical state 
according to their needs. It is clear by now that this 
mechanical state is not simply predefined by biochemical 
processes, but that cells can be characterized as complex 
materials with nontrivial nonlinear mechanical feedback. 
Cell mechanics also couples back to the physiological re- 
sponse of cells [2J, for example to the spreading behavior 
[3] or even to stem cell differentiation [3]. In this con- 
tribution, we aim at elucidating the fundamental physics 
providing cells with their unique material properties. We 
make use of the fact that the mechanical response of the 
cell can be ascribed to the cytoskeleton [5], a complex 
biopolymer network. Within the constituents of the cy- 
toskeleton the dominant response is contributed by the 
actin cortex, a semiflexible polymer network connected 
by weak, reversible crosslinks. It is the microstructure of 
this actin cortex that we use as starting point for mod- 
eling. More precisely, we base our discussion on a model 
for transiently crosslinked biopolymer networks, called 
the inelastic glassy wormlike chain (inelastic Gwlc) [6]. 
From the inelastic Gwlc, which is numerically still quite 
demanding, we extract a reduced constitutive model that 
is more tractable both analytically and numerically. We 
restrict ourselves to reversible ("no-slip") on-off kinetics, 
but find the model nevertheless capable of describing a 
large part of the known experimental data. For illustra- 



tion, we evaluate the response to a standard deformation 
protocol. 



II. THE INELASTIC GWLC 

The inelastic glassy wormlike chain (inelastic Gwlc) 
|6J is a phenomenological mean-field type model that de- 
scribes a test polymer fluctuating against a background 
network. The slowing-down of the conformational dy- 
namics of the test polymer by the surrounding network 
is phenomenologically represented by a stretching of its 
mode spectrum. More precisely, the single mode relax- 
ation time r n of the nth eigenmode of a free polymer is 
modified according to 



r n exp [£(A„/A - 1)] 
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In the spirit of a mean-field model, £ and A are below 
interpreted as the characteristic free energy barrier of, 
and the average contour distance between, weak tran- 
sient bonds, respectively, and A„ is the wavelength of the 
nth mode of the polymers's transverse contour fluctu- 
ations rx(s,t). Inserting the modified relaxation spec- 
trum, Eq. ([T|), into the expression for the mechanical re- 
sponse of an isolated wormlike chain in solution yields 
a highly successful phenomenological parametrization of 
the equilibrium response of sticky biopolymer networks 
[7] and even of living cells jS]. For example, the sus- 
ceptibility of a test polymer of bending rigidity k to a 
transverse point force is given by [9] 
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where /a is the Euler buckling force. 
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We found the above interpretation of the mathematical 
expression in terms of physical network parameters par- 
ticularly useful for non-equilibrium processes, in which 
the bond network is driven out of equilibrium so that 
the average bond distance A deviates from its force-free 
reference value Aq according to 
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(3) 



where v quantifies the fraction of broken bonds and v$ 
is the respective equilibrium value. For simplicity, we 
model v by a first-order kinetic equation with two possi- 
ble states, bound and unbound, characterized by on and 
off rates fc + and fc_ , respectively [5] 

u(t, f) = - {k+(f) + M/)} "(*) + M/)- (4) 

The rates are prescribed by Kramers' theory [TO] and an 
exponential force dependence according to Bell's model 
QI], i.e. 



and 



k + (f) = ke u ~ A ^f , 



fc_(/) = ke Ax " f , 



(3) 



(G) 



with k ~ Tq e~ £ . They depend on the height Ek^T 
of the the energy barrier separating bound and unbound 
state, the relative binding affinity e , and the widths Axt 
and Ax u of the bound and unbound state, respectively. 
In the following, the thermal energy k^T is set to one for 
convenience. 



III. SIMPLIFIED CONSTITUTIVE EQUATIONS 

The crucial feature of the inelastic Gwlc is that the 
mutual interaction of the "stiff" polymer response and 
the "soft" bond network generates a huge variety of phe- 
nomena, depending on the particular stimulus. Here, we 
are interested in the long-time quasi-plastic response, and 
not in short-time effects and the particular shape of the 
relaxation spectrum, which we discussed elsewhere [3]. 
We therefore refrain from reproducing all the details of 
the full model and concentrate on those aspects that are 
most important in this respect, namely a nearly expo- 
nential strain stiffening and the softening due to the re- 
versible breaking of weak bonds. This leads to a reduced 
simplified formulation of the inelastic Gwlc that allows 
for the analytical derivation of constitutive equations for 
transiently crosslinked biopolymer networks. 

To arrive at these equations, we cast the model equa- 
tions into a form that emphasizes common aspects of 
plasticity theory. As a starting point, the system size L 
is related to the initial size Lq by a scaling function a nl 
that characterizes the nonlinear compliance in terms of 
the dimensionless bond fraction v, 



The inelastic Gwlc model predicts a monotonic decrease 
of the material stiffness with the bond fraction [B], i.e. the 
more bonds are broken, the more susceptible the material 
becomes to deformations. For the sake of the argument, 
we approximate the precise functional form by a simple 
reciprocal dependence of the nonlinear susceptibility on 
the bond fraction, 



[/,!/(/)] = ^<K/>0) 



(8) 



Here, we have introduced the (nonlinearly) elastic com- 
ponent 6. It only depends on the equilibrium bond frac- 
tion uq under zero external force, while the dependence 
on the current bond fraction v is isolated into the pref- 
actor. Defining the inelastic length C as the part of the 
extension that is not due to elastic deformations, 

C = L i/ /i/, (9) 

the actual system length can now be written as 

L = C-5. (10) 

Stressing the similarities to plasticity theory, we can 
associate a "rest force" T with the inelastic length, de- 
fined by the condition C = 0. Given our knowledge of 
the bond kinetics, the rest force is easily calculated. Ac- 
cording to Eq. (bl, £ = z> = 0. The rest force T is 
therefore determined by the following condition 



F ■ v(t)\jr = 0. 

By Eqs. Q-Q, for a given force /, v = 
(1 + e (^»+^)f-uy\ yielding 



(Ax b + Ax u )T = U + In 



(11) 
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(12) 



where we solved for f and replaced v using relation QQj) . 
In the context of bond kinetics, T is the force that would 
have to be maintained to render the current bond fraction 
stationary. Taking a time derivate, the rest force is seen 
to inherit its dynamics from the inelastic deformation 
rate £, a phenomenon commonly denoted as kinematic 
hardening. However, in contrast to the common linear 
kinematic hardening (T cx £), bond breaking leads to 
the phenomenology of logarithmic kinematic hardening, 
T oc C/C. Note that T is finite for C > LqVq and that, 
by definition, J- = for C = L a . 

The corresponding flow rule can be found by decom- 
posing the total mechanical force / = T + Af into the 
rest force T and an overstress A/, and inserting it into 
the force-dependent bond kinetics, Eqs. Q-^. Elim- 
inating T using Eq. (12 1, and expressing the result in 



terms of the inelastic length C, the flow rule is given by 
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where we introduced the abbreviation e = Ax\,/(Axb + 
Ax u ). At any instant t, Eq. (TL31) uniquely determines 
C, given the overstress history — A/(i). Note that the 
equations derived above do not depend on the particular 
choice of the elastic response 5. 

To summarize this section, we found a set of constitu- 
tive equations from an approximate model based on the 
molecular structure of the cytoskeleton, describing the 
inelastic mechanical response of a transiently crosslinked 
biopolymer network. We can now proceed to evaluate 
the equations and to numerically examine the responses 
to a simple deformation protocol. 



IV. QUASI-PLASTIC RESPONSE 

We consider the response of the model to a deforma- 
tion ramp T = L/Lq = To = const., followed by a 
plateau T = and an inverse ramp, t = —To (Fig. [T] 
lower panel, inset). The elastic contribution of the inelas- 
tic Gwlc is well approximated by exponential elasticity, 
f(S) cx (5 — l^e 7 '^ -1 ) , which has been observed for a 
multitude of biomaterials [T5] . For simplicity, we also ne- 
glect the viscous component of the viscoelastic polymer 
response, as it is not essential for the following argumen- 
tation. The qualitative effects presented in the following 
do not depend on these technically motivated simplifica- 
tions. 

The force-displacement curves exhibit rate-dependent 
hysteresis (Fig. [I] upper panel), which in our model is 
a signature of dissipation caused by inelastic reversible 
bond breaking. In the limit of an infinitely slow deforma- 
tion, the bonds are always in equilibrium and the hystere- 
sis vanishes (solid lines in the upper panel of Fig. [I}, since 
we dismissed the viscoelastic hysteresis present in the 
Gwlc model and real biopolymer materials, here. For 
low to moderate rates (the time scale is set by the intrin- 
sic time scale T$e £ of the bonds), all force-displacement 
curves share characteristic features. Most prominently, 
the response to a linear ramp is characterized by two 
approximately linear regimes, emphasized by the dashed 
lines in the upper panel of Fig. [T] The first linear regime 
can be interpreted as an elastic response. The second 
regime is explained by a quasi-plastic deformation with 
slowly increasing rest force a signature of the appar- 
ent kinematic hardening identified above [see Eq. ( 12 1], 



0.05 



where we found that the rest force depends logarithmi- 
cally on the inelastic deformation, suggesting the notion 
of "logarithmic kinematic hardening". Corresponding de- 
viations from the linear force-displacement curve become 
discernible for deformations with amplitudes much larger 
than 10% (not shown). The phenomenology is reminis- 
cent of experimental results obtained for living cells |13j . 

The interpretation of an elastic and an inelastic regime 
is substantiated by comparing the time-dependent rest 
force T to the total mechanical force / (Fig. [I] lower 
panel). While during a ramp, the total mechanical force 
initially strongly diverges from the rest force, it quickly 
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Figure 1. Response of the simplified schematic model for 
biopolymer networks, Eqs. ( 10 1- ( 13 1 , derived from the in- 
elastic Gwlc, to a piece- wise linear deformation protocol 
(loading- waiting-unloading, see inset of lower panel). The 
relative deformation amplitude is AF = 10%, and ramp 
and plateau durations are AF/Tq. Upper panel: force- 
displacement curves [deformation rates FoToe £ = 0.004 (solid 
lines), 0.64 (dotted lines), and 1.6 (dot-dashed lines)]; dashed 
lines indicating elastic and apparent kinematic hardening 
regimes, respectively, reminiscent of what is observed for liv- 
ing cells |13| . Lower panel: time-dependent force for the high- 
est deformation rate (dot-dashed lines) compared to the time- 
dependent rest force T (solid lines); forces are normalized by 
the natural force scale Ax^ (in units of ksT) of the bound 
state. 



settles on a course parallel to the rest force, consistent 
with a constant elastic contribution. The elastic con- 
tribution is relaxed upon halting the deformation. The 
force characterizing the transition between elastic and 
inelastic regime can be interpreted as an effective yield 
threshold. For infinitely slow deformations, the total me- 
chanical force equals the rest force, and no predominantly 
elastic regime is present. In other words, also the ef- 
fective yield threshold depends on the deformation rate. 
This observation sets our biopolymer network apart from 
usual models for hard solids. The reason for this behav- 
ior is that, by construction, the bonds will always yield if 
the stimulus is slow compared to the bond-opening time 
scale. Only for sufficiently fast deformation, an initial 
predominantly elastic response can be obtained. 
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V. CONCLUSIONS 



Starting from a minimal model of a transiently 
crosslinked biopolymer network, we derived quasi-plastic 
constitutive equations exhibiting logarithmic kinematic 
hardening. In contrast to "truly" plastic materials, the in- 
elastic deformations are rooted in the inelastic, reversible 
softening due to the transient breaking of weak bonds, 
as presumed by the cell rheological model of the inelastic 
Gwlc. 

From the perspective of biomechanics, the present 
work may lead to an intuitive understanding of the me- 
chanical properties of biopolymer materials and cells. It 
also provides a simple but accurate model to simulate 
the large-scale behavior of the materials, e.g. using finite 
element methods. From the perspective of materials sci- 
ence, our work sheds light on a new class of materials, 
which can bear a kind of strain that is both recoverable 
and dissipative, as opposed to the "usual" reversible elas- 
tic and the irreversible plastic strain. This recoverable 



inelastic strain bears resemblance to the quasi-plastic- 
elastic (QPE) model recently proposed in relation with 
DP steel alloys [2]. In contrast to the QPE strain, how- 
ever, our recoverable inelastic strain is rate-dependent, 
due to the underlying slow bond dynamics. It is an in- 
triguing question whether the QPE strain might emerge 
from the bond-breaking approach in some special limit- 
ing case. As an outlook, we would like to mention that 
it would be straightforward to extend our model to ac- 
count for true plastic strains [13] (by associating some slip 
with each bond breaking event) as well as for the phys- 
iologically important internally generated active stresses 

USUI]. 
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